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I. INTRODUCTION 

The complexity of the hadronic mass spectra induces the possibihty that existing and newly observed hadrons 
can possibly be interpreted as molecular states (or hadronic molecules). Such an interpretation is possible, when 
the mass of the hadronic molecule rriH lies slirfitlybelow the threshold of the corresponding hadronic pair H1H2: 
niH < ruHi + (for review see e.g. Refs. In the light meson sector, possible candidates for hadronic 

molecules are the scalar mesons ao(980) and ^(980) treated as KK bound states [^. Including the heavy flavor 
meson sector other possible molecular states can arise. For example, the scalar and axial charm _D*n(2317), -Dsi (2 460) 
and bottom S*o(5725) and Ssi(5778) mesons can be treated as DK, D*K, BK and B*K bound states 0, S iM El, 113 , 
respectively. Other candidates for a hadronic molecule interpretation are the X(3872) as a D^D*^ + charge conjugate 
(c.c) bound state, r(4260) as a DDi - c.c. and ^(4415) as a DlDso{2?,n) + c.c. bound state In the baryonic 

sector, the most popular candidate for a hadronic molecule is the negative-parity l/2~ resonance A(1405) considered 
as a NK bound state [s|. Also, there are candidates in the heavy baryon sector, e.g. the charmed baryon Ac(2940)+ 
recently discovered by the BABAR Collaboration 14] which can be treated as a Z?*°p bound state flSj . 

In the current manuscript we focus on the scalar charm-strange meson Z?*q(2317), which was discovered just a 
few years ago by the BABAR Collaboration at SLAC in the inclusive Dfir^ invariant mass distribution of e+e" 
annihilation data [l^. The nearby state Dsi{2A60) with a mass of 2.4589 GeV decaying into D*tt° was observed 
by the CLEO Collaboration at CESR 17]. Both of these states have been confirmed by the Belle Collaboration at 
KEKB |18i]. From interpretation of these experiments it was suggested that the £'*g(2317) and Z3si(2460) mesons 
are the P-wave charm-strange quark states with spin-parity quantum numbers ~ 0+ and = 1+ states, 
respectively. In the following the Belle ^1^ and the BABAR [20] Collaborations observed the production of £'*q(2317) 
and -Dsi(2460) in nonleptonic two-body B decays together with their subsequent strong and radiative transitions. 
Taking into account existing experimental information on the properties of D*q{2317) and £'si(2460) mesons [2lj, one 
can conclude that the respective = 0"*" and = 1"*" quantum numbers are now established with high confidence. 

The next important question concerns the possible structure of the £'*o(2317) and Dsi(2460) mesons. The simplest 
interpretation of these states is that they are the missing jg = 1/2 (the angular momentum of the s-quark) members 
of the cs L — I multiplet. However, this standard quark model scenario is in disagreement with experimental 
observation since the Z3*q(2317) and Dsi(2460) states are narrower and their masses are lower when compared to 
theoretical (see e.g. discussion in Ref. [sj). Therefore, in addition to the standard quark-antiquark picture alternative 
interpretation of the Z3*q(2317) and I?si(2460) mesons have been suggested: four-quark states, mixing of two- and 
four-quark states, two-diquark states and two-meson molecular states. Up to now different properties of the D*g(2317) 
and Dsi{2A6Q) mesons (masses, strong, radiative and weak decay constants and widths) have been calculated using 
different approaches [3]j[l3]"[i3,[23"lS3- quark models, effective Lagrangian approaches, QCD sum rules, lattice 
QCD, etc. 

In present paper we will consider the strong D*q ^ Ds + tt'^ and radiative D*q D* -I- 7 decays of the £'*q(2317) 
meson using an effective Lagrangian approach. The approach is based on the hypothesis that the D*q is a strong bound 
state of D and K mesons. In other words we investigate the position that Z?*g meson is a (DK) hadronic molecule. 
The coupling of the D*q meson to the constituents {D and K mesons) is described by the effective Lagrangian. The 
corresponding coupling constant go'gDK is determined by the compositeness condition Z — js^ . [ssl . [Sa ]. which 
implies that the rcnormalization constant of the hadron wave function is set equal to zero. Note, that this condition 
was originally applied to the study of the deuteron as bound state of proton and neutron [s^ . Then it was extensively 
used in the low-energy hadron phenomenology as the master equation for the treatment of mesons and baryons as 
bound states of light and heavy constituent quarks (see Refs. [H, [s^, [H, [13] )■ In addition this condition was 
used in Ref. 61] in the application to glueballs as bound states of gluons. Recently the compositeness condition 
was used to study the light scalar mesons qq and /o as KK molecules A new impact of the DK molecular 
structure of the D*q(2317) meson is that the presence of u{d) quarks in the D and K meson gives rise to a direct 
strong isospin-violating transition D*q DsTt^ in addition to the decay induced hy rj — n'^ mixing considered before 
in the literature. We show that the direct transition dominates over the 77 — tt" mixing transitions. The obtained 
results for the partial decay widths are consistent with previous calculations. By analogy one can treat the second 
charm narrow resonance Dsi(2460) as a {D*K) molecule and the possible corresponding bottom counterparts - the 
states B*q(5725) and Bsi(5778) - as BK and B*K bound states, respectively. The calculation of the properties of the 
Z?si(2460), -6*0(5725) and i3si(5778) mesons goes beyond the scope of the present paper and we relegate this issue to 
a forthcoming paper. Also in near future we plan to consider two-body B-meson decays and semileptonic processes 
involving D*o(2317) and 1531(2460) in the final state. 

In the present manuscript we proceed as follows. First, in Section II we discuss the basic notions of our approach. 
We derive the effective mesonic Lagrangian for the treatment of charm and bottom mesons £'*o(2317), Dsi(2460), 
B*q(5725) and i?si(5778) as DK, D*K, BK and B* K bound states, respectively. We discuss how to determine 
the corresponding coupling constant between the hadronic molecule and its constituents using the compositeness 
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condition. In Section III we consider the matrix elements (Feynman diagrams) describing the strong and radiative 
decays of the £'*q(2317). We indicate our numerical results and discuss various limits, such as the local case and the 
heavy quark limit. In Section IV we present a short summary of our results. 



II. APPROACH 



Molecular structure of the D*q*(2317) meson 



In this section we derive the formalism for the study of the D*^{2317) meson as a hadronic molecule - a bound 
state of D and K mesons. First of all we specify the quantum numbers of the Z?*g^(2317) mesons. We use the current 
results for the quantum numbers of isospin, spin and parity: I{J^) — 0(0+) and mass — 2.3173 GeV [IJ. Our 

framework is based on an effective interaction Lagrangian describing the coupling between the 15*0(2317) meson and 
their constituents - D and K mesons: 



^D'Jx) = g^, D:^{x) dy<i>D'Jy')D'ix + WKy)K{x-WDy) + H.c. (1) 



The doublets of D and K mesons are defined as 

D-[^l), ^-('^o), (2) 

the symbol T refers to the transpose of the doublet D. In particular, the assumed molecular structure of D^q and 
D*Q states is: 

\D:+) = \D+K°) + \D'K+) , \D:^) ^ \D-K°) + \D°K-) . (3) 

The correlation function $_d*q characterizes the finite size of the D*q(2317) meson as a (DK) bound state and depends 
on the relative Jacobi coordinate y with x being the center of mass (CM) coordinate. Note, the local limit corresponds 
to the substitution of ^d'^ by the Dirac delta- function: $£)*^(j/^) 6'^{y). The kinematical variables wd and wk 
are defined by 

iriD rriK , 

WD = ■ , WK = ■ , (4) 

1710+ mK mD+mK 

where mo and ttik are the masses of D and K mesons. The Fourier transform of the correlation function reads 

^o'Jy') - l^,--'''^o:,i-p'). (5) 

Any choice for <I>_D*g is appropriate as long as it falls off sufficiently fast in the ultraviolet region of Euclidean space 
to render the Feynman diagrams ultraviolet finite. We employ the Gaussian form 

^d:M)^^M-pI/^1'J, (6) 

for the vertex function, where pe is the Euclidean Jacobi momentum. Here ^d^„ is a size parameter, which 
parametrizes the distribution of D and K mesons inside the D*q molecule. 

The D*^DK coupling constant go'^^ is determined by the compositeness condition [s^, HEH^i which implies that 
the renormalization constant of the hadron wave function is set equal to zero: 

= l-S'c.^/m^.J=0, (7) 
where S^. {'m-'o- ) = 5^ , n'^^* {'m-'o* ) is the derivative of the D^^ meson mass operator described by the diagram in 

As we already stressed in Introduction, this condition was originally applied to the study of the deuteron as a 
bound state of proton and neutron [54| . Then it was extensively used in low-energy hadron phenomenology as the 
master equation for the treatment of mesons and baryons as bound states of light and heavy constituent quarks 
(see Refs. [H, [13, [H, |H, [60]). In Ref. [6l| this condition was used in the consideration of glueballs as bound 
states of gluons. Recently the compositeness condition was applied to the study of the light scalar mesons aq and 
/o as KK molecules [o*]. To clarify the physical meaning of this condition, we first want to remind the reader 
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that the renormahzation constant Z^', can also be interpreted as the matrix element between the physical and the 
corresponding bare state. For -^d'q = it then follows that the physical state does not contain the bare one and is 
solely described as a bound state. The interaction Lagrangian Eq. ([T]) and the corresponding free parts describe both 
the constituents {D and K mesons) and the hadronic molecule {D*q), which is taken to be the bound state of the 
constituents. As a result of the interaction the physical particle is dressed, i.e. its mass and its wave function have to 
be renormalized. The condition Zd*^ = also guarantees that there is no double counting for the physical observable 
under consideration: the D*^ meson interacts with other hadrons and gauge bosons only via its constituents. In 
particular, the compositeness condition excludes the direct interaction of the dressed charged particle (like D*^ 
mesons) with the electromagnetic field. Taking into account both the tree-level diagram and the diagrams with the 
self-energy and counter-term insertions into the external legs (that is the tree- level diagram times (Zd'^ — 1)) one 
obtains a common factor Zn*^ which is equal to zero j55„ ^ST,, ^] . 



B. Effective Lagrangian for strong and radiative decays of D*q*(2317) 

Now we turn to the discussion of the lowest-order diagrams which contribute to the matrix elements of the strong 
isospin-violating decay D*^^ — > DsTt'^ and the radiative decay D*^ — s- D*j. To the strong decay two types of diagrams 
contribute: the so-called "direct" diagrams of Fig. 2 with 7r°-meson emission from the D^*^K meson loops and the 
"indirect" diagrams of Fig. 3 where a tt" meson is produced via ry — tt" mixing. Note, that the second mechanism based 
on ?7 — 7r° mixing was mainly considered before in the literature. Originally, it was initiated by the analysis based on 
the use of chiral Lagrangians [13, [H, HI] where the leading-order, tree-level D^qDsTt^ coupling can be generated only 
by virtual 77-meson emission. During the last years different approaches have been applied to the D*q — > DsTt^ decay 
properties using the — tt'^ mixing mechanism. In our approach the D*q meson is considered as a DK bound state 
and, therefore, we have an additional mechanism for generating the D*f^Ds'K^ transition due to the direct coupling 
of and iff*) mesons to 7r°. In particular, in the isospin limit (when the masses of the virtual and K'^*^ 
mesons in the loops are degenerate, respectively) the pairs of diagrams related to Fig. 2(a), 2(b) and Fig. 2(c) and 2(d) 
compensate each other. Only the use of physical masses for the and K^*^ mesons gives a nontrivial contribution 
to the — > Ds7r° couphng of order 0[5), where 

5 - m^(., ± - 77i^(.) - TO^(.) ± - m^(., (8) 

is the parameter of isospin breaking. Therefore, the contribution of the diagrams of Fig. 2 is of the same order as the 
one related to Fig. 3 involving rj — tt^ mixing, where the 77 — tt" transition coupling (filled black circle) is counted as 
0{S). 

The diagrams contributing to the radiative decay D*^^ D*^^ are shown in Fig. 4. The diagrams of Figs. 4(a) 
and 4(b) are generated by the direct coupling of the charged and K'^ mesons to the electromagnetic field after 
gauging of the free Lagrangians related to these mesons. The diagrams of Figs. 4(c) and 4(d) (so-called contact 
diagrams) are generated after gauging of nonlocal strong Lagrangian ([1]) describing the coupling of D*q mesons to its 
constituents - D and K mesons. The diagrams of Figs. 4(e) and 4(f) arise after gauging the strong D*DK interaction 
Lagrangian containing derivatives acting on the pseudoscalar fields. Finally, the diagrams of Figs. 4(g) and 4(h) 
describe the sub-process where the D*^ converts into the D* via a DK loop followed by the interaction of the D* with 
the electromagnetic field. Note that an analogous diagram where the D* meson interacts with the electromagnetic 
field and then converts into the D* vanishes due to the transversity condition for the on-shell vector meson D*, i.e. 
P/j (p) = 0. Details of how to generate the effective couplings of the involved mesons to the electromagnetic field 
will be discussed later. 

After the preliminary discussion of the relevant diagrams, now we are in the position to write down the full effective 
Lagrangian Cctt for the study of strong D*q — )■ DsTt'^ and radiative D*q —> D*j decay properties. For convenience we 
split >Ccfr into an isospin-symmetric part >Cinv and an isospin-symmetry breaking part >Cbrcak: 

Ccff{x) = Cinv{x) + Chiciik{x) , (9) 

where C-mv is given by a sum of free meson parts £froc and the interaction parts C-mt- 

Cinvix) = C[rcc{x) + Cint{x) ■ (10) 

We use the standard free meson Lagrangian involving states with quantum numbers = 0^, 0^ and 1^: 



^frco(2^) — ^ ^froo(^) J 



(11) 



5 



where 



Cf,,,ix) = -D: + ix){D + mi,,jD:~{x), (12) 

- DHx)in + ml)D{x)~Dt{x){n + mljD-{x), (13) 

£L,{x) = K;\x){g'^-'[a + ml,] - d^dnKix) + D;\x){g'^''[n + ml,] - d'^dnOtix) 

+ D: + {x){g^^'^[D + ml,] ~ 9^9")^?: -(x) . (14) 

Here □ = 9^9^, tt is the triplet of pions, Df and D*^ are the pseudoscalar and vector charm-strange mesons, 
respectively. The doublets of vector mesons D* and K* are given by 

D* =( nZ) ^ K*=( fll ] . (15) 



^ D*+ J ' V ^ / 

In our convention the isospin-symmetric meson masses of the iso-multiplets are identified with the masses of the 
charged partners [2lj |: 

m^ = m^± = 139.57018 MeV , = mK± = 493.677 MeV , mK* = mn-. ± = 891.66 MeV , (16) 

mo = m^± = 1.8693 GeV, mo- = tod.± = 2.010 GeV. 

The masses of the iso-singlet states are j2l| : 

m,, = 547.51 MeV , mo, = m^,* = 1.9682 GeV , (17) 

mo* — m^*± = 2.112 GeV, "^d'^ = mjj,± ~ 2.3173 GeV. 

The interaction term C{at{x) will be discussed later. First we would like to write down the isospin-breaking term 
'Cbroak, which includcs the mass corrections of the neutral mesons containing u or d quarks and the r] — tt'^ mass 
mixing [g^, HI]: 



>Cbrcak(x) = SC^'ix) + SC^'ix) + L,^{x) , (18) 



where 



SC^'ix) = ^-^[ti\x)Y ^5kK''{x)K\x)+5dD\x)D''(x), (19) 
5C^{x) = -5k' K*J'{x)K*°^{x) - 5d' Dl''{x)D*'^^'{x) , (20) 
C,^{x) = B!!M^!^^0(x)r,(x), (21) 

where m„ and are the u and d current quark masses, B is the condensate parameter. Here 5m are the isospin- 
breaking parameters which are fixed by the difference of masses squared of the charged and neutral members of the 
iso-multiplets as: 

5m = rnli± - m\jo , tomo = mjQ-o . (22) 

The set of rriMO is taken from data [2l| with: 

m^o = 134.9766 MeV , Wj^o = 497.648 MeV , m^.o = 896.0 MeV , (23) 
m^jo = 1.8645 GeV, TO£,.o = 2.0067 GeV . 

Eqs. p2| - (fT4l) . (|19p and ((20| define the free meson propagators for scalar (pseudoscalar) fields 

iDM{x-y) = {Q]TM{x)M\ym ^ 1 e-'=(--^) ^,,(fc) , (24) 
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where 



'M 



le 



and vector fields 



where 



{2it)H 



-'''(^-y) (fc) 



(25) 



(26) 



(27) 



In the following calculations it will be convenient to expand the propagators of the neutral mesons £'"(13°), K'^\K''^), 
]j*0(^£)*o^ and K*'^{K*'^) in powers of the corresponding isospin-breaking parameters as: 



DMoik) = 



1-5 



d 



M- 



1-5 



d 



DM±{k) + 0{5li), 



DZ^Ak) + 0{5l,,: 



(28) 



The interaction Lagrangian includes the strong and electromagnetic parts 

Ant(.T)-/:rnI(-T^)+^rntW, 



(29) 



as already apparent from the previous discussion related to Figs. 2-4. The relevant strong part of the effective La- 
grangian contains the following terms: the Lagrangian Cd^^ ([I]) describing the coupling of the _D*q meson to its 
constituents and yPP-type Langrangians, describing the interaction of vector mesons with two pseudoscalars: 

+ Cd*d,k{x) + Ck*d,d{x) + Cdidk{x) . (30) 
Let us specify the VPP interaction Lagrangians occurring in Eqs. ((30|) . In general they can be defined as: 

J^VP^pM = <?v^p,p, V^,{x)Pi{x)iT P2{x) + H.c. (31) 
To be consistent with the definitions occurring in literature, we use the following form of the particular Lagrangians: 



C-D*D-q{x) 
jO-K'KTrix) 
C-K*K-n{x) 
CwD^Kix) 

J^k*d,d{x) 

Cd'Dk{x) 



-^^^^D*Ux)TT:ix)ilj^ Dix) + H.c. 
2%/2 ^ 

_^^.t(^)^(^),5^^(^) +H.C. 

^^^^K*^{x)TTT{x)id^ K{x) + H.C. 

^-J^K;;{x)rj{x)td'K{x) +H.C. 

g^.^^^D*J{x)K{x)iT D-{x) + H.c. 
g,,,^,Kf{x)D{x)iTD-{x) + H.c. 
9^,,^D:-{x)D'^{x)iT K{x) + H.c. 



(32) 

(33) 

(34) 

(35) 

(36) 
(37) 
(38) 



where summation over isospin indices is understood and A d B = AdB — BdA. 

The couplings (?o*D,r ^^'^ 9k'k^ fixed by data for the strong decay widths D* Dtt and K* Ktt. In 
particular, the strong two-body decay widths r{D*+ £)°7r+) and r(A'*+ K°Tr+) are related to g^,^^ (ei. [ssj 
and g^.^^ as 



T{D*+ D"t:^ 



9\ 

Dtt 



247rm^^+ 



pi 



9\ 



(39) 
(40) 
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where P^^ is the three-momentum of 7r+ in the rest frame of the decaying vector meson V . Using data for the 
corresponding strong decay widths one deduces: g^y,^^ = 17.9 [11] and 5^.^^ — 4.61 [21i] . 

The coupUng constants 90*07,^ obtained in the context of heavy hadron chiral perturbation theory 

(HHChPT) [66| . The couphngs g^,^^ and g^.^^ are expressed (and then related) in terms of a universal strong 
coupling constant g involving heavy (vector and pseudoscalar) and Goldstone mesons and in terms of the leptonic 
decay constants Fp: 

mo. /TT mw [2 

= -^5V2, 5^.^^ = -^5y-, (41) 

where F^^ = 92.4 MeV and — 1.3 F^^. From Eq. (|^T|) and using g^,^^ — 17.9 we deduce the value of gj^.^^ with 

5.... = = 7.95. (42) 

The coupling constant g^,^^ can be related to g^.j^^ using the unitary symmetry relation: 

5...-,, = ^5..., =6.14. (43) 
Again, as in the case of g^,^ ^ j, we include in couplings the relation to the corresponding decay constants F-j^ and 

Fr,. 

The couphng constants 5^.^ ^ and ff^.^^ have been estimated using the QCD sum rule technique in Refs. [67l.[68t. 
These couplings are important for the evaluation of the dissociation cross section of J/^ to kaons (see, e.g. discussion 
in Refs. [i^, Here we use the results of Ref. [13]: S^.d k = 2. 02 and g^,„i^ — 1.84. The couphng 5^.^ „ can 

also be related to 5^.^ ^, using SU{A) symmetry arguments: g^^.j^ „ ~ 9d*d.i< ~ 2.02. 

The relevant electromagnetic part has three main terms: 

a-(^) = C^Z^'\x) + C^Z^'\x) + CZ^'\x) . (44) 
The first term describes the local coupling of charged D-, K- and _D* mesons to the electromagnetic field 

Cf'^(^) = ieA^{x)i^D-{x)T D+{x)+K-{x)T K+{x) 

- d:-^x) Td*+{x) + \d:-"{x) dcD*+'^{x) + ii?r^(.T) Td:+{x)^ . (45) 

The term is generated after gauging of the free meson Lagrangians using minimal substitution: 

Qf^M^ (a^ zp iev4^) Af=^ . (46) 

The terms ^I'nt^^^ ^'^'^ '^in?*''^^ generated due to the gauging of the strong Lagrangians ([5T|) and (P) containing 
derivatives acting on the charged fields. Note, that the correlation function describing the nonlocal D*qDK cou- 

pling, is a function of and, therefore, both Lagrangians pip and ([1]) are not gauge-invariant under electromagnetic 
C/em(l) transformations and should be modified accordingly. 

To get the second term we replace all derivatives acting on the charged fields by the covariant ones using minimal 
substitution (as is the case for gauging the free Lagrangians). The term in relevant for our calculation in 

contains the coupling of the vector D* meson to D, K and the photon field with 

')(x) = eg,,,,A^ix)D:;{x)[D°ix)K+ix) D+{x)K°{x)] + H.c. + ••■ (47) 

The g aug in g of the nonlocal Lagrangian of Eq. ^ proceeds in a way suggested in Ref. [7l[ and extensively used in 
Refs. |57l l58l|. In particular, to guarantee local invariance of the strong interaction Lagrangian, in Cf^l each charged 
constituent meson field (i.e. and meson fields) are multiplied by the gauge field exponential resulting in 



^str+cm(3)(^^ = ^sV (^) j dy <PDt,{y') [e-^^'^^+^-y'^'''^ D+ {x + WKy)K° {x - wnv) 

+ D°{x + WKy)e'"''^''~"''"y^'^'^^K+{x-WDy)\ + B..C. (48) 
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where 

/(x,2/,P) - / dz^A^^iz). (49) 



V 



For the derivative of the path integral (UHl) we use the path-independent prescription suggested in Refs. [7l| 

d 

Urn dx''—I(x,y,P) ^ \im \I(x + dx,y,P')-I(x,y,P)], (50) 

dx^^—^0 OX^ dxt'-^O 

where path P' is obtained from P when shifting the end-point x by dx. Use of the definition (jSO]) leads to the key 
rule 



-^I{x,y,P) = A^{x), (51) 



which in turn states that the derivative of the path integral I{x, y, P) does not depend on the path P originally used 
in the definition. The non-minimal substitution (|48l) is therefore completely equivalent to the minimal prescription. 

In the calculation of the amplitudes of the radiative D*q D*^ decay, in Eq. (|48p we only need to keep terms 
linear in A^, that is the four-particle coupling D*qDK^. Hence, the third term contributing to the electromagnetic 
interaction Lagrangian is given by 

x+wkV 

^Z^'\^) = -^eg^,^^D:,-{x) Jdy^n-Jy')!^ J dz^A^^{z)D+ {x + WKy)K\x - woy) (52) 

X 

x—woy 

dZf,A''{z)D°{x + WKy)K^{x - woy) \ + H.c. H 



Concluding the discussion of the effective interaction Lagrangian we stress that all couplings occurring in the diagrams 
contributing to the decays D*q DsTt'^ and D*q — > D*j are explicitly fixed, except g^, discussed in the following. 



C. Analysis of the DIqDK coupling g^. 



Finally, we discuss the numerical value of the model-dependent constant g^, . In terms of a general functional form 
of the correlation function '^di^ the coupling constant g^, is given by: 

C30 oo 



where 



gl, (47rAr,* )2 J J (l + ai+a2)3 

"sO 



2 2 ^^^D*„ 

z = ^.j^ai + Ai^a2 - — ■ , (54) 

1 + ai 4" Q!2 

R = aia2 -I- aiWjj + a2Wj^ , = -. • 

An* 

One should stress that coupling constant g^, remains finite when we remove the cutoff A^i*^^ oo (local limit). 
A finite result is obtained, because the derivative of the D*q mass operator is convergent, i.e. the loop integral is 
J d'^k/k^ when the correlation function ^D*g is removed (or equal to one) at A^,.^ oo. However, in the calculation 
of transition diagrams (like in Figs. 2-4) we deal with divergent integral and, therefore, we need the correlation function 
to perform the regularization of the occurring loop integrals. Now the question is how sensitive our results are to a 
variation of ^d'^- First, we look at the coupling constant g^, . In the limit ^d'^ ^ oo it is given by 



. In^ - 1 + ^ 7^ " ^ y arctan^ (55) 
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where z± = m'jjt ± {mj~, — itt-'k) and 

A = A(m|,.^, to|,, TO^) = m%,^ + m% + m% - 2m\,.m\, - 2m\,.m\ - 2m\,m\ (56) 
is the Kallen function. 

For checking purposes we also analyze the coupling g^, in the heavy quark limit (HQL), where the masses of D 

and -D*Q mesons together with the charm quark mass rric go to infinity. In the HQL the D meson in the D*q molecule 
fixes the center-of-mass, surrounded by a light K meson in analogy to the heavy-light Qq mesons. For the nonlocal 
case the result for in the HQL is: 



1 I f da 



Ol, (47rmc)2 J 1 + /i|a 



K'M)- (57) 



The HQL result for the local case is: 



9 

mi 



In^ . (58) 

(47rmc) mj 



Now we compare our numerical results for the coupling constant g^, in different regimes: 1) nonlocal case (NC); 
2) local case (LC); 3) nonlocal case + HQL (NCHQL) and 4) local case + HQL (LCHQL). When we deal with the 
nonlocal case we proceed with the Gaussian correlation function $£)»^(z) = exp(— z) and vary the scale parameter 
^D*g from 1 to 2 GeV. For the charm quark mass we use the averaged result of the PDG ^2l\: rric = 1.25 GeV. For a 
convenience we attach a corresponding superscript to g^, to indicate the specific regime. 

We get the following results: the coupling g^^ varies from 11.26 GeV at A^,.^ = 1 GeV to 9.90 GeV at A^,*^, = 2 GeV. 
The coupling g^^ is expressed only in terms of physical meson masses with the resuh .g^^ = 8.98 GeV. The couphng 
^NCHQL yj^j.-gg f^Q^-^ '^Q 22 GeV at A^,.^ = 1 GeV to 11.52 GeV at Ad;^^ = 2 GeV. Finally, we have 5^?"*^^ = 11.52 
GeV. All results for g^, are quite close to each other with a typical value for g^, of about 10 GeV which is consistent 
with preceding calculations done in other theoretical approaches. In Table 1 we compare our result for the D*.^DK 
coupling constant to predictions of other theoretical approaches (we use a compilation of the results done in Ref. [67^). 



D. Extension to other possible hadronic molecules 

We end this section with a comment concerning the extension of the derived framework to the study of other 
hadronic molecules. This can be done in a straightforward fashion. The starting point is the construction of an 
effective Lagrangian describing hadronic molecules as bound states of its constituents. In particular, for the case of 
the charm-strange meson I?si(2460) and for the possible partners in the bottom sector (5725) and i3si(5778) the 
simplest Lagrangians have the form: 

^D., {x) = g^^^ D;/{x) I dy <^>n^, (y2) D*/{x + w^^, y) K(x - w^.^y) + H.c. , (59) 
Cb'Jx) = g,,^ B*J{x) J dyfB'Jy^) B\x + w^^y) K{x - w,^y) + H.c. , (60) 

^bAx) = 9s^,B°r{x) jdy<^BAy^)Bl\x + w^^,y)K{x~w^.^^y) +11.C., (61) 

where wij = mi/{ini + rrij), gM and $jv/ are the coupling constants (fixed from the compositeness condition) and 
correlation functions. The doublets of _B^*) and _B(*)t mesons are defined as 

= ( ) : B^*'^ t ^ - 0) _ (g2) 

The molecular structure of Df^, B*§ , 13* § , and B°^ is: 

\D+) ^ \D*+K") + \D*^K+) , \D-^) = \D*-K") + \D*^K-) , 

\B*^) ^ \B+K-) + IB^K"^) , IB^O) = \B-K+) + {B'^K'^) (63) 

|B,"i> = \B*+K-) + \B*"K") , = \B*-K+) + \B*"K°) . 

The calculation of decay properties of -Dsi(2460), B*q{5725) and _B5i(5778) mesons goes beyond the scope of the 
present paper and we relegate this issue to a forthcoming paper. 
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III. STRONG Dto -> Dsn° AND RADIATIVE D*o Dty DECAYS 

In this section we discuss the numerical results for the D*q Dgir^ and -D*o — > D*^ decay properties. As we 
already stressed in the preceding section two types of diagrams contribute to the amplitude of the strong decay 
D*Q Dg-K^: the "direct" diagrams of Fig. 2 and the "77 — 7r° mixing" diagrams of Fig. 3. The "direct" diagrams occur 
due to the DK molecular structure of the D*g meson, while in the two-quark picture they are forbidden according to 
the Okubo-Zweig-Iizuka rule. In the framework of our approach this is not the case, since D and K mesons contain 
nonstrange quarks. The total contribution of the "direct" diagrams starts at order 0(5), where 5 of Eq. ([8|) is the 
generic parameter of isospin breaking. Hence, the "direct diagrams" are of the same order in the isospin-breaking 
counting scheme as the "77— 7r° mixing" diagrams, and, therefore, both types of diagrams should be included. To clarify 
this mechanism we present our results our results for two cases: 1) "full calculation" (Full) and 2) "leading-order" 
(LO), i.e. restricting to first order in isospin-breaking 0{5). 

It is convenient to write the matrix element describing the _D*q — > DsTt'^ transition as a sum of the added contribu- 
tions of the diagrams in Figs. 2 and 3: 

M{D*o ^ DsTT°) = MdniD*o ^ Dsn") + Mn^Dto ^ ^«^°) (64) 

with 

MdiriDla -> i^.^°) = 9Di„D,.o , (65) 

M,^UD:o - D,n") = , (66) 

nis ~ m 4 

where m = (m,i + md)/2 and (m^ — mu)/{ms — rh) — 1/43.7 (see e.g. Ref. t62j). In the derivation of Eq. (j66p we use 
the masses of tt" and rj meson in leading order of the chiral expansion. The total effective D*i^Dsn" coupling, denoted 
as GD'gD^n includes both contributions of the set of diagrams of Figs. 2 and 3 with 

GD'gD.TT = G^'JoD.Tr + GS'oD.TT : (67) 

sO all aO s sa I _ ^ 4 

In terms of Gd'^d^i) the D*q DsTt" decay width reads as: 

no:, ^ D,n) = p;o , (69) 



Stttti^, 







where P*o = A"'^/^(to|,,^ , m|)^ , m^o)/(2mD«^) is the three- momentum of the decay products. 

The matrix element describing the D*q — > I?* 7 transition can be written in the manifestly gauge- invariant form 

M^,{D:„ ^ D:^) - e GD;,Di^ {g^up'q - p'^q.) , (70) 

where p' and q are the D* and photon four- momenta and p — p' + q is the D*q momentum. Here G 0*^0*1 is the 
effective D*qD*^ coupling constant and the _D*q D*j decay width is given by 

r{D:, -> d:j) = a Gi.^^ot, p;' (71) 

where 
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(72) 



is the three-momentum of the decay products. 

Now we present the numerical results. First we discuss the contributions of the different diagrams of Figs. 2 and 3 
to the effective coupling Gdi^d^-k- With a typical value for the scale parameter of ^di^ = 1 GeV we get the following. 

In the Full calculation, 

Gd-„d.. = 146.6 MeV, G^J^^^, = 104.5 MeV , Ggf^^,. = 42.1 MeV , (73) 
Ggr^^, = 40.9 MeV , G^'j;^fJ^ = 63.6 MeV , G^i^'^^f^ = 7.9 MeV , G^^^^,^*^ = 34.1 MeV . 



11 



and in the LO calculation, 

Gd*„d,^ - 145.4 MeV , G^'J^d,^ = 103.4 MeV , G^i^^^^ = 42.0 MeV , (74) 
Gt'foZ'^ = 40.2 MeV , G^'^.^^fJ^ = 63.2 MeV , G^f'^^f - 7.9 MeV , G^f '^^^*^ = 34.1 MeV , 

where the superscripts DD*K and KK*D relate to the diagrams of Figs. 2(a), 2(b), 3(a), and 3(b) and 2(c), 2(d), 
3(c), and 3(d), respectively. The direct diagrams dominate over the mixing diagrams by about a factor of 2. The 
results for the decay width (total result and partial contributions of the different diagrams) are as follows. 

In the Full calculation, 

r(i5:„^I?,^)=46.7KeV, TiDl^ 
T{Dla DsTif"^^^''^ = 3.6 KeV, 
r(D*o -> D.tt)™'^'-"^'^^' = 0.1 KcV, 

In the LO calculation, 

r{D% DsTt) = 46.6 KeV , r(D,*o ^ Dsirf" = 23.6 KeV , r(D,*o ^ I?s7r)""'^ = 3.9 KeV , (76) 
r{D*o ^ Dsirf"'^^''^ = 3.6 KeV , r{D*o Dsirf"''^'^''" = 8.8 KeV , 
r{D;^ -> Li.Tr)'"'''-^^'^ = 0.1 KeV , T{Dlo ^ L'.tt)""'^'-^^*'" = 2.6 KeV . 

From Eqs. ([75)1 and ([75)1 it is evident that the restriction to the leading-order in isospin breaking is very good 
approximation to the full calculation (both sets of results practically coincide with each other). We would like to 
stress that the strong decay width T{D*q Dgir) is enhanced in a molecular picture for the D*^ meson as compared 
to the quarkonium interpretation due to the inclusion of the direct tt*^ coupling to the DD* or KK* meson pairs. 
This enhancement is particularly present, since the "direct" mode dominates over the "mixing" mode. 

On the other hand, when turning to the heavy quark limit the contribution of the "direct" mode becomes much 
smaller, about 0.4 KeV, while the one of the "mixing" decreases less to about 1.4 KeV. The total result for the decay 
width is an order of magnitude smaller as in the full dynamical case with T{DIq — > DgTi) — 3.3 KeV. From the results 
obtained in the HQL we make the following conclusions: first, in the HQL the "mixing" mode dominates over the 
"direct" mode. This result is consistent with heavy hadron ChPT by conception (restriction to the "mixing" mode) 
and numerically (the result for the width is of the order of a few KeV). Second, we have a clear explanation why 
in the HQL the "direct" mode is suppressed. The reason is that the isospin breaking effects due the difference of 
heavy mesons occurring in the loop are of next-to-leading order in the 1/toc expansion, i.e. they are of the form 
(5£)(.) / (mcAij*^). Numerically these factors are not so small when compared to the isospin-breaking factors /A|)* 

arising from the mass differences of kaons if^*-'. We conclude from our results, that the heavy quark limit is not a good 
approximation for the isospin- violating strong decay D*q — > Dgir^ since some of the important isospin-breaking effects 
are missing. In addition, taking in general the HQL in the charm sector is not necessarily a good approximation 
because of the relatively small mass of the charm quark. In contrast we show below that for the radiative decay 
D*o Z?*7 the HQL works well. 

In Table 2 we present our results for the decay width T{D*q — > Dgi:) including a variation of the scale parameter 
Ai3*^ from 1 to 2 GeV (increase of A/j*^ leads to an increase of the width) and compare them to previous theoretical 
predictions. 

Now we turn to the discussion of the radiative decay D*^ D*^. By construction, using a gauge-invariant 
and Lorentz-covariant effective Lagrangian, the full amplitude for this process is gauge- invariant, while the separate 
contributions of the different diagrams of Fig. 4 are not. It is important to stress that the diagrams of Fig. 4 fall 
into two separately gauge-invariant sets: one set includes the diagrams of Figs. 4(a), 4(c), 4(e), and 4(g) (with loops 
containing virtual and mesons), generated by the coupling of Z?*q to the D+ and constituents. Another 
set contains the diagrams of Figs. 4(b), 4(d), 4(f), and 4(h) (with loops containing virtual D° and mesons) with 
the coupling of to £)" and . 

For convenience wc split each individual diagram into a gauge-invariant piece and a reminder, which is non- invariant. 
One can prove that the sum of the non-invariant terms vanishes due to gauge invariance. In the following discussion 
of the numerical results we will deal only with the gauge-invariant contribution of the separate diagrams of Fig. 4. 
Another important feature of the D*g D*^ amplitude is that the effective coupling Go'^D^-y is ultraviolet (UV) 
finite. In the Appendix we discuss the local limit that is when we remove the cutoff with A^.^ — > oo in the correlation 



DgTrf'' = 23.7 KeV , r(D*o ^ i^s^r)™ 3.8 KeV , (75) 
r(L>:o ^ DsT^f"^'^'^'^ = 8.8 KeV, 
r(i:i:o ^ DgT:)'^^''^'^'^ = 2.5 KeV. 
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function ^d*„- Again, the separate contributions of the diagrams of Fig. 4 to Gdi^D'-i contain divergences which 
compensate each other. In the Appendix we discuss this issue in detail. 

First, we show the results for the effective coupling constant Gd'^d^j- the total result and partial contributions of 
the different diagrams of Fig. 4 (marked by 4(a), 4(b), etc.). In the analysis of the electromagnetic decay £)*q — » 0*^ 
we restrict to the isospin limit, i.e. we do not include the isospin-breaking effects in the meson masses and proceed 
with the masses of the charged particles. In the isospin limit the diagrams of Fig. 4(e) and 4(f) compensate each other 
(and therefore do not contribute to the total amplitude), while the diagrams of Fig. 4(g) and 4(h) are equal to each 
other. For a value of A^-^ = 1 GeV we get 

Gd'^d-i = 0-093 GeV"^ , Gj^.^jj,^ = -0.030 GeV"^ , Gd-„D'^^ = 0-089 GeV"^ , (77) 
^D^oi^ti = 10"' G^V-i , G^i^j,.^ = 0.002 GeV-^ , G^f.^^,^ ^ G%^J,,^ = 0.016 GeV"! . 
The corresponding results for the decay width D*q — > Z)*7 are: 

r(£>:o^D:7) = 0.47KeV, (78) 

r(D:o ^ D:7)4- = 0.05 KeV , r(D:o ^ D:7)4b = 0.43 KeV , 

r(D*o ^ 0*^)^" 6 X 10"^ KeV , r{Dl„ Dl^f)'^'^ = 2 x lO"" KeV , 

r(D:o i?:7)4g ^ r(z?:o ^ d:^)^^ ^ 0.02 Kev. 

From the results it is clear that the contact diagrams of Fig. 4(c) and 4(d) are strongly suppressed, these diagrams 
are kept to guarantee gauge invariance. The main contribution comes from the diagram of Fig. 4(b) where the photon 
couples to the K^. The diagram of Fig. 4(a) is relatively suppressed as ^ [niK /mof' . 

The sum of all the diagrams is ultraviolet finite and, therefore, the cutoff parameter can be removed with A^)*^ oo. 
In the local approximation (LC case) for the radiative decay width we get the following results for the coupling 
constant Gd* D"y and the decay width r(Z)*o [Here we only deal with the gauge-invariant parts of diagrams 

of Figs.4(a),4(b), 4(g), and 4(h).]: 

Gd'^^d^j = 0.110 GeV-i , 

Gd' D". = -0.038 GeV"^ , G% r,,^ = 0.093 GeV"^ , 

^ Gj,^^J,,^ = 0.055 GeV-i , (79) 

and 

T{D:o^Dtj)=OM KeV, 

T{D*o D*'y)'^^ = 0.08 KeV , T{D*„ 0*^)*"° = 0.47 KeV , 

r{D;^ ^ D*j)^^ = r{D:^ ^ d:^)^^ ^ 0.04 KeV . (80) 

The LC results are larger than for the nonlocal case (NC case) choosing A^*^ = 1 GeV. 

Finally, we consider the HQL to this process. In the NCHQL case the diagrams of Fig. 4 relatively scale as: 



-(4a . ,-(46 . ^4c . r<4d . ^ig{h) _ 1 . i . ^ . 1 . 

^ D* r)*'Y • *-^D* r)*'v • *-^D* r)*'v * ^n* d*'v * ^n* n*'v — * ^ * — o • ^ • 

'^so'^sl ^so'^sJ ^.o'^sJ J-'sO^e-y J^sO^s^ rUn TO^ TOc 



Therefore, the leading order contribution arises from the diagrams of Fig. 4(b) and 4(d), resulting in 

Gd'^„D'^^ = 0.114 GeV-i , 

Gf;,^o,^ = 0.053 GeV"^ , G')^,^o,^ = 0.061 GeV"^ , (82) 
and the corresponding results for the decay width of 

r(D:o ^ ^?:7) = 0.71 KeV, 

r(L>;o ^ Dl'yf'^ = 0.15 KeV , r(£)*o ^ Dl-if"^ = 0.20 KeV . (83) 
Finally, in the LCHQL case the diagrams of Fig. 4 relatively scale as: 

^D'^oDti ■ ^D'^oDti • ^dI^-t ^ : 1 : 1 (84) 
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Therefore, the leading order contribution comes from the diagram of Fig. 4(b) with 

— — In — : 

and 



Gd,„d,, = In^ = 0.160 GeV- , (85) 



r(i^:o-> 7^:7) = 1.41 KeV, (86) 

where the coupUng g^, is given by Eq. ([58|) . 

In Table 3 we summarize our results for r(Z?*Q D*j) for all four cases (NC, LC, NCHQL and LCHQL) including 
a variation of the scale parameter Ad»^ from 1 to 2 GeV (an increase of ^d*„ leads to a larger value for the width). We 
also compare to predictions of other theoretical approaches. Our results have a negligible dependence on the parameter 
^D*f^ and are also in good agreement with previous calculations. Also, within a factor of 2 our results for the different 
considered cases are in good agreement. Hence for the radiative decay D*q — > Z?*7 the local approximation (LC) and 
HQL are reasonable approximations. 



IV. SUMMARY 



We studied the new charm-strange meson £)*q(2317) in the hadronic molecule interpretation, considering a bound 
state of D and K mesons. Using an effective Lagrangian approach we calculated the strong D*q DgTi^ and radiative 
D*Q iD*7 decays. A new impact of the DK molecular structure of the D*o(2317) meson is that the presence of 
u{d) quarks in the D and K meson loops gives rise to a direct strong isospin-violating transition D*q — > D^tt^ in 
addition to the decay mechanism induced by r/ — 7r° mixing as was considered before in the literature. We showed 
that the direct transition dominates over the 77 — tt" mixing transition. Our results for the partial decay widths are 
summarized as follows: 



T{D*Q D,tt) ^ 79.3 ± 32.6 KeV ["Full" calculation] , 

r(D*o ^ D^tt) = 79.6 ± 33.0 KeV ["LO" calculation] , 

(87) 

r(D*o -> D*j) = 0.55 ± 0.08 KeV ["NC" case] , 

r{Dlg DIj) = 0.66 KeV ["LC" case] , 

r{D'*() 13*7) = 0.94 ± 0.23 KeV ["NCHQL" case] , 

r{D'*() 13*7) = 1.41 KeV ["LCHQL" case] . 

The ratio R = r(L'*Q D*j/r{DlQ Dsn) ~ 10^^ satisfies the current experimental upper limit of i? < 0.059 [2l|. 

For the case of the strong decay the application of the heavy quark limit (HQL) gives a significant suppression of 
the direct mode. The contributions of the isospin-breaking effects associated with the mass-difference of Z3^*^ mesons 
have an extra factor AD*^/mc and, therefore, are formally of higher-order in the l/rric expansion in comparison to 
the isospin-breaking effects associated with the mass difference of K^*") mesons. However, numerically the factor 
Ad*(,/™c is of order 1, leading to the result that the HQL is not a suitable approximation for the isospin-violating 
decay I3*g ^ DsTt°. 

In the case of the radiative decay D*q D*^ we have another situation and the different limiting cases (local 
limit, heavy quark limit) considered give more or less a similar description of the physical quantities GD'^D^-y and 
r(-DIo ~^ ^*sl) (see the results of Table 3). Here our conclusion is that in the context of a molecular interpretation 
the decay width r(D*o D^j) is of order 1 KeV as was previously predicted before by other theoretical approaches. 
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APPENDIX A: MATRIX ELEMENT OF THE RADIATIVE DECAY Dlo D^j 

Here we discuss the matrix element of the radiative decay D*q ^ D*^ in the local approximation (when the cutoff 

in the D*Q meson correlation function is removed with ^d'„ ^ oo) and for the nonlocal case. 

As we mentioned before the on-shell matrix element describing the D*q D*"f transition can be written in the 
manifestly gauge-invariant form 

M^,u{Dl„ ^ Dli) = e Gdi.D'^^ {a^..p'q - p'^^q.) ■ (Al) 

In the local approximation the following diagrams of Fig. 4 contribute to this matrix element: diagrams of Figs. 4(a), 
4(b), and 4(e)-4(h). As we stressed in Sec. Ill, two sets of all the diagrams are separately gauge-invariant: the set 
of Figs. 4(a), 4(e), and 4(g) and set of diagrams related to Figs. 4(b), 4(f), and 4(h). For illustration we take one set 
[Figs. 4(a), 4(e), and 4(g)] and prove gauge-invarianco by using dimensional regularization (DR) for the separation of 
the divergent pieces which finally cancel each other. 

The structure integrals (we drop the occurring coupling constants) corresponding to the diagrams of Figs. 4(a), 4(e), 
and 4(g) are given by: 



In diagram Fig. 4(a) 



In diagram Fig. 4(e) 



In diagram Fig. 4(g) 



J (27r)^i [m% -{k+ p)2] [ml -{k+ p')'] [ml - k^] ' ^ ' 

T^t = -9^.. j (2^^ [^2^ _ (fc + p)2] [^2^ _ fc2] (A3) 



where 



^ -g-^+pVMl,. r rPk {2k + p')p 

ml, - p2 J (27r)^z [ml - {k + p)^] [ml - {k + p')'] [ml - ' ^ ' 



r^.a = - {P + pX + ^{P + P')u +^{P + P% . (A5) 



Next using the Feynman a-parametrization and the master formula of DR 

jDi, ( u2\M 



I 



d^'k i-k'r ^ 1 T{D/2 + M)r{N-M-D/2) d/2+m-n 
(27r)^i [A - (47r)^/2 r{D/2)T{N) 



we get: 



T^'. = ^ j da{l-2a)\n\nK + ^^Jml ^d^-P'^^- P^-) j da{l - 2a)\nl^DK + 0{D - A) , 



(A6) 



1 



1 3 

- -^{9^.p'q-p',q.) j d:'a5(l ~ ^ + 0{D - 4) , (A7) 

'^-^ 

1 

= -|^{4^+^"^" + ^'w} + i^ / da\n^.nK+0{D-A), (A8) 



(A9) 
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where 

^DDK = A3{mD,mK) = mjy{l - as) +ml^a3 ~ m^jj.^aias ~ rnjj.a2a3, (AlO) 
Adk = A2{mD,mK) — m%{l — a) + mj^a — mj2i.^^a{l — a) . 

From Eqs. (|A7P - (|A9p one can see that in the sum of the diagrams of Figs. 4(a), 4(e), and 4(g) all divergences and 
non-gauge invariant pieces cancel each other. Taking I? ^ 4 we write down the final result of: 

1 3 1 

T',:^'"^'' = ^i9,^P'Q-p',Q^)[-J d'as(^l~J2^^) + ^ / da{l-2a)\nAnK^ (All) 

By analogy we prove the gauge invariance for the sum of the diagrams of Figs. 4(b), 4(f), and 4(h): 

1 3 1 

^M^'^^'' = ^i9,.p'q~p',q.)[J d^as(^l-Y,a?j - ^ J da{l ~ 2a)\nAKD^ (A12) 



'-^ " 



where Akkd = A^^mfcrni)) and Akd — A2(mK,rnu)- It is easy to show that the second terms in Eqs. (jAll[) and 
(|A12[) are equal to each other by changing the variable a to 1 — a. Therefore, the total result for the effective coupling 
constant Gd* D*-y in the local case is: 



3 1 



G.,,.,, = [jSasil - E a) aiaaj^ - ^} + ^ j Ml - 2.)lnA..} . (A13) 

In the nonlocal case the gauge invariance can be proved based on a method developed e.g. in Ref. [5^. For this 
purpose in particular we split the contribution of each diagram into a part which is gauge invariant and one which is 
not: we use the following representation for the four-vectors with open Lorentz indices /i and v. 

P' -P±;,+l'^, (A14) 
p = p^.p, +P 

such that p'^.gqfj, = and P±.pip'i, = 0. Expressions for diagrams containing only ±- values are gauge invariant 
separately. It is easy to show tnat the remaining terms, which are not gauge invariant, cancel each other in total. 
Note that this method works perfectly both for on-shell and off-shell amplitudes. 
The coupling constant Go'^D'-y in the nonlocal case is given by 

Go'^^Dij = '-^g^2 ^D^^Dij , (A15) 

where iD'gD^-y is the structure integral containing the contributions of the diagrams in Figs.4(a)-4(d), 4(g), and 4(h): 

i—a.b.c,d,g,h 

oo oo oo 

r4a _ ^ J J J ^"^^"^^"^ 





^5*o^.*7 = -J^J J J (l + ai23)^ («i+^^)(«3 + u>i^) hd<l>'o.JznDK)], 



oo oo oo 



^^'o^.-T = i/// ti + Ztr i-d^'D'J^KKD)], (A16) 





1 





1 oo oo 

Aj 4 2 /" I f f daida2 , n r ji/ / m 

'■Dl^Ol'-^ = J^'^D J "-^^ J J (-1_^Q^^)4 [WDa2-WKai) [-d^D*„[ZKD)\, 
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where 



oo oo 

rig _ Tih ^ f f daida2 , \ ^ t \ 



=° 



ai23 = ai + a2 + as , ai2 — ai + a2 , 

ZdDK = Zz{lJiD, IJ-k) , ZkKD = Z3{llK, I^d) , ZdK = Z2{lJ-D , IJ-k) , ZkD = Z2{lJ-K, IJ-d) , 

zsifJ-i, 1^-2) = /uf ai2 + Maas + mI,. «;iW2 - '^^ '^'^^ (ai + ^2) + lJ,%,a2) , (A17) 

su -|- Cti23 * * 

/ \ 2 2 / 2 2 /-I \\ C^2 ~t~ / 2 2 / /-i \ \\ 

22(/^i,M2j = Ml"! +/^2"2 + (m/j. t + -^jj'^^iW'a - -r- (/x^,. w;2i + M£>*(w'2(l - ij , 



,2 „, , ,,2 „, , ,,2 /„_„.,_. (ai +t«if)(a2 +W'Z3)^ WM 



-ZP = + Mifa2 + HD'r. I WDWK 7- , MM 
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Table 1. Coupling constant g^, 
The range of values for our results is due 
to the variation of An» from 1 to 2 GeV. 



Approach 


9o:^nK (GeV) 


Ref. [51] 


2.5 - 3.8 


Ref. [72] 


5.068 


Ref. [73] 


5.5 ± 1.8 


Ref. [67] 




Ref. [38] 


6.0 - 7.8 


Ref. [49] 


9.3t|I 


Ref. [74] 


< 9.86 


Ref. [n] 


10.203 


Our results; 




NC case 


9.90 - 11.26 


LC case 


8.98 


NCHQL case 


11.52 - 16.22 


LCHQL case 


11.52 



Table 2. Decay width of DIq ^ Dsn°. 
The range of values for our results is due 
to the variation of Adj^^ from 1 to 2 GeV. 



Approach 


TiDto ^ Ds7T°) (KeV) 


Ref. [43] 


6 ± 2 


Ref. [27] 


7 ± 1 


Ref. [22] 


10 


Ref. [31] 


16 


Ref. [28] 


21.5 


Ref. [45] 


32 


Ref. [42] 


39 ± 5 


Ref. [32] 


15 - 70 


Ref. [26] 


10 - 100 


Ref. [33] 


129 ± 43 (109 ± 16) 


Our results: 




Full case 


46.7 - 111.9 


LO case 


46.6 - 112.6 



Table 3. Decay width of D^q —f Dl^ . 

The range of values for our results is due 
to the variation of An» from 1 to 2 GeV. 



Approach 


V{D% ^ D't'y) (KcV) 


Ref. [31] 


0.2 


Ref. [27] 


0.85 ± 0.05 


Ref. [40] 


1 


Ref. [44] 


1.1 


Ref. [50] 


1.3 - 9.9 


Ref. [33] 


< 1-4 


Ref. [28] 


1.74 


Ref. [22] 


1.9 


Ref. [3^ 


4-6 


Ref. [32] 


21 


Our results: 




NC case 


0.47 - 0.63 


LC case 


0.66 


NCHQL case 


0.71 - 1.17 


LCHQL case 


1.41 



